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Abstract 



A general stability condition for plasma-vacuum systems with re- 
sistive walls is derived by using the Frieman Rotenberg lagrangian 
stability formulation [Rev. Mod. Phys. 32, 898 (I960)]. It is shown 
that the Lyapunov stability limit for external modes does not depend 
upon the gyroscopic term but upon the sign of the perturbed poten- 
tial energy only. In the absence of dissipation in the plasma such as 
viscosity, it is expected that the flow cannot stabilize the system. 

PACS: 52.30.Jb, 52.35.Py, 02.90.+p 
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1 Introduction 



In a previous theorem [T] a general instability statement was proved con- 
cerning static plasmas surrounded by resistive walls. This statement is re- 
lated to the problem of resistive wall modes (RWM) in the magnetohydrody- 
namic (MHD)literature. In particular, the impact of plasma rotation on the 
mode has been investigated theoretically in Refs.[2]-[8] and experimentally in 
Refs.[9]-[TT]. See also Ref. [12] for mode calculations in a cylindrical case and 
the recent review paper of Ref. [13]. In the present paper we show that the 
general Lyapunov analysis can be extended to stationary plasmas surrounded 
by resistive walls even if their resistivity is time-dependent as investigated in 
Ref.[H]. 

In section II we discuss the Lagrangian stability investigation of Ref . [15] 
in the context of free boundary displacements through an interface between 
plasma and vacuum. The Lyapunov functional is constructed in section III 
while conclusions and discussions are left to section IV. 

2 Stability equation 

The linearized Lagrangian displacement vector £ about an ideal MHD equi- 
librium with flow obeys the equations derived in Ref. [15]: 

Po 'i + 2p v • V£ + F£ = 0, (1) 

where po is the equilibrium mass density and vo the flow velocity. The anti- 
symmetric 2p vo • V operator is sometime called the "gyroscopic" operator 
and F is given by 

F£ = V(7PoV-^-Vpo-B-Q)+B-VQ+Q-VB+V-(po€vo-Vvo-pov vo-V^), 

(2) 

where po is the equilibrium pressure and B the equilibrium magnetic field. 
Q = V x (£ x B) is the perturbed magnetic field. The operator F of Eq.(2) 
contains in particular flow terms so generalizing the corresponding operator 
of the static case [TB] . 

After fixing the gauge by annihilating the scalar potential, the vector 
potential A in the region outside the plasma is governed by Ohm's law 

V x V x A = -o-A (3) 
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where a is the space and time dependent conductivity in the outer region. a 
vanishes for a true vacuum. 

At the interface between the plasma and the outer region, we must fulfil 
the condition 

nxA=-(n-£)B (4) 

where B is the equilibrium magnetic field at the vacuum side of the interface 
whose normal vector is n. A has to vanish at infinity or be normal to 
superconducting walls, if there are any. Equations (1) to (4) reduce to the 
equations (1) to (3) of Ref.[T4"] if the flow vanishes. 

The gyroscopic operator 2poVo- V is antisymmetric and the scalar product 
2 / £ • (A) v o " V)£dr as a volume integral vanishes for n • £ = at the plasma 
boundary. If n ■ £ ^ at the interface the scalar product vanishes also if the 
mass density vanishes at the boundary. To make sure that this is true we 
use 

V x (£ x p v ) = Pov • V£ - £ ■ V(p v ) - p v V • £, (5) 

V(£ ■ p v ) = p v • V£ + p v x V x £ + £ ■ Vp v + £ x V x p v . (6) 

Note that the term £V • (poVo) has vanishing contribution to Eq.(5) because 
of the equilibrium continuity equation. Adding equations (5) and (6) we 
obtain 

V(£-p v ) = 2poVo-V£-poVoV-£-Vx(£xpoVo)+pov xVx£+£xVxpoVo. 

(7) 

After taking the scalar product of Eq.(7) with 77 and using 

V • [77 x (£ x p v )] = (£ x p vo) - Vxt7-?7-Vx(£x p v ), (8) 

V • [(77 • p v )£] = (rj ■ p v )V • £ + £ ■ V(rj ■ p v ), (9) 

we integrate over the plasma volume and transform the volume integrals of 
the divergences to surface integrals at the interface to obtain 

J r) ■ (Pov ■ V)£dr = - J £ ■ (p v • V)r)dr (10) 

if either n • £ = n ■ 77 = for internal modes or poVo = at the interface for 
external modes. 
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3 Lyapunov functional 



Taking the scalar product of Eq.(l) with £ and Eq.(3) with A and integrating 
over the plasma volume and the vacuum region respectively we obtain 

d(K + W) , ■ 

— — -(ctA, A) (11) 

where 

(e,Fa + (A,VxVxA) or 



w 



K 



{f>ok,k)p 



the subscripts p and or being related to the plasma and the outer region 
respectively. The parentheses denote the scalar products 



Jp 



vn,i*)or = / vn,v>F 



or 



Since a is positive, the quantity K + W has a negative time derivative for 
external modes according to Eq.(ll). This allows us to apply the method of 
Lyapunov, which states that if K + W has no definite sign, then the system is 
"Lyapunov unstable". For internal modes with A = the condition is only 
sufficient. The Hermitian form W could be minimized in the same way as the 
5W of the energy principle [16] for static MHD equilibria, but no numerical 
code for minimizing W has been developed yet. For a straight tokamak with 
homogeneous axial flow and without walls W is indefinite if it is indefinite 
in the static case. This follows from the fact that the stability property 
does not depend upon the inertial frame chosen. Therefore, according to 
the above Lyapunov study the introduction of resistive walls cannot stabilize 
the plasma-wall system with finite constant plasma velocity. One should 
note, however, that Lyapunov stability is equivalent to spectral stability if 
all modes of the spectrum are considered. This is not in agreement of the 
result of Ref. [Bj. In general, the entire spectrum of the plasma- wall system 
is not as easy to investigate as the determination of the sign of W. 

Eq.(ll) reduces to Eq.(4) of Ref. [TJ [H] if the velocity terms in F vanish. 
It extends the Lyapunov functional introduced in Ref. [TJ [TJ] to the present 
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general flow case. This means that if W is positive definite the system is 
stable, but if it is indefinite the system is unstable. So the gyroscopic operator 
does not affect stability of external modes but growth rates, eigenvalues and 
eigenfunctions only. 

4 Discussion and conclusions 

The question of stability of external modes for toroidal ideal plasmas with 
flows surrounded by resistive walls is decided by the sign of the symmetric 
potential operator only, which is a drastic simplification so that numerical 
investigation of stability becomes only slightly more difficult than in the 
static case for which many codes have been developed. 

Moreover the Lyapunov stability property holds also for internal modes 
if a dissipation operator is added to the gyroscopic operator in (1) since this 
addition would not change the sign of the right hand side of (11). However, if 
we add a viscosity to the physical system under investigation then (1) would 
be modified by such an operator acting on £ but also by another operator act- 
ing on £ which represents the "circulatory" forces as considered in Ref.jlT]. 
The operator on £ is proportional to the viscosity itself but the "circulatory" 
operator goes with viscosity times the velocity vo. Only the latter operator 
can modify the stability property, so a minimum velocity is needed for this 
operator to be able to stabilize the system. However, the quantitative inves- 
tigation of such situations is difficult as explained in Ref.[T7] and Ref . [18] . 
We should also mention that the plasma viscosity is a tensor known only in 
the collisional regime while a convincing derivation for the weak collisionality 
regime of fusion plasmas has not been done yet. Therefore, it may turn out 
that the stability of hot plasmas must be treated by kinetic equations like 
Vlasov or Fokker-Planck for which the present stability procedure does not 
apply. 

Finally, if for physical interpretations of the experiment, eigenmodes and 
growth rates of (1) to (4) (and possibly nonlinear effects) are needed, then 
the investigation of such situations will be much more difficult especially for 
toroidal plasmas. Numerical techniques for the search of complex eigenvalues 
must be developed in the first place. (See e.g. [19]). 
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